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In the bearingless induction motor (BIM), the rotor of BIM is supported by 
the magnetic force, the radial position control technique is studied in this 
paper, an approximation and simplification are used to create a linearized 
model of suspension system. The objective of this paper is to design a control 
system that overcomes the limitations of the traditional controller and to 


improve the control performance and achieve high stability, even if it is 
exposed to external disturbance and internal disturbance (changes in load or 
speed). In order to solve this problem, the proposed controller is based on the 
fuzzy-I controller that combines the positive features of both fuzzy logic and 
the integration. In comparison with the conventional PID, the results show 
that the fuzzy-I controller reduces the peak-to-peak rotor deviation by 23% 
under effect of external disturbance force and 31% under effect of speed 
variation and 9% under effect of load variation. 
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1. INTRODUCTION 

A mechanical bearings motor is unable to meeting the requirements for high speed and long-term 
operation. As a result, the magnetic bearing motor was created. However, the magnetic bearing motor has 
many disadvantages such as complex structure, high energy consumption for the magnetic suspension, low 
efficiency, difficulty of over-speed and so on [1], [2]. The bearingless motor is the best solution for this problem 
where the number of the inverters, dimensions of the motor, and the associated costs are reduced [3], [4]. 

There are many types of the bearingless motor, the bearingless induction motor is one of the most 
important of these types. The BIM is a hybrid of conventional motor and magnetic bearing, The BIM consists 
of two set of windings are imbedded in stator slots, the electromagnetic torque windings to generates rotation 
torque and the suspension control windings to control the suspension force [5]—[8]. A magnetic levitation force 
is generated by the interaction of the air-gap magnetic flux generated by two sets of the stator windings, used 
to control the axial suspension force of the rotor [1]. Compared with magnetic bearing, the bearingless motor 
has many advantages, like higher critical speed, shorter shaft, smaller in size, and lower cost [9]. In a 
bearingless induction motors, the volt-Ampere required for axial levitation force control is very small compared 
with the volt-Ampere of torque windings [4]. This result indicates possible cost reduction in the bearingless 
induction motor and demonstrates it is less expensive. The induction motor has many advantages over other 
types of motors such as its robustness, low cost and its simplicity [10]. 

The important problems and requirements for controlling BIM is to generate an electromagnetic force 
in order to prevent the collision between the stator and the rotor that makes the rotor rotates in the center of 
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stator bore. A sliding mode controller with variable structures of BIM was proposed by [11] to implement the 
dynamic decoupling control with excellent performance and enhance the stability of the system. Direct torque 
and the suspension forces control for BIM using active disturbance rejection controller was proposed [12], the 
simulation and testing findings demonstrate that the ADRC technique has excellent robustness, can effectively 
reduce the influence of the load disturbance and also reduce overshoot. Direct levitation force control method 
using model prediction has been proposed by [13], the modeling and experiment results demonstrate that the 
suggested control technique can successfully increase the suspension stability of the rotor also enhancing the 
system's ability to resist disturbances. The non-singular fast terminals with sliding mode control was proposed 
in [14], The experiment and simulation result demonstrate that the proposed controller tracks the values of the 
radial displacement and speed quickly, also, the rotor is reaching to the steady state position immediately and 
achieving stability. 

The goal of this paper is to design control system to overcome the traditional controller limitations 
and to improve control performance and achieve high system stability even when the motor shaft is subjected 
to external disturbances. The proposed control system in this work is based on the fuzzy-I controller: this 
controller combines the features of both fuzzy and integration. fuzzy logic is a nonlinear control system, it 
doesn’t require exact mathematical models and the integration helps to reduce the steady state error. 


2. MATHMATICAL MODELS OF BEARINGLESS INDUCTION MOTOR 

The BIM contains two groups of stator winding are embedded in the same slots. In this work, the 
following structure is considered: four poles winding pı to generate the electromagnetic torque and two poles 
winding p2 to generate the rotor suspension forces as shown in Figure 1. The two groups must achieve two 
important conditions in order to generate stable suspension force: @1=@2 and P;=P2+1, where œ and @2 
represent the angular frequencies of torque and suspension windings respectively. 

Figure 2 Shows that the flux linkage in air-gab Y, and ‘¥2 are generated by the current of torque 
windings la and suspension windings 2a respectively, ‘¥; and ‘¥2 are in same direction on the left side of air- 
gap, ¥ı and Y, are in opposite direction in right side of air-gap, as a result the magnetic flux density in left side 
is greater than the magnetic flux density on the right side, then the radial levitation force F is generated in the 
negative X axis. If the currents in radial suspension windings are reversed, a radial levitation force is generated 
in the positive X axis. In the same way the radial levitation force in Y direction is controlled by supplying 
current in the suspension windings 2b which are perpendicular to the suspension windings 2a [15]. 


Roto 
Torque winding 


Rotor shaft 
Suspension winding 
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Figure 1. Structure of BIM 
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Figure 2. Principle of radial the suspension forces generation 
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2.1. Electromagnetic torque and speed of BIM 

The magnetic flux of the axial suspension windings and torque windings induce a current in the rotor, 
however, the current induced by the suspension windings and torque that caused by it can be ignored because 
it is very small with respect to the current induced by torque windings, so the magnetic flux of the air-gap can 
be written as [14]-[16]: 


Wid =Lim Chas + har) (1) 
Vig = Lim (igs + lig) (2) 


where: Yiq and Yiq are the flux linkage in air-gab in d-q coordinates system respectively, Lim is mutual 
inductance for torque windings, itar, iigr, iias and ijgs are the rotor and stator current components for torque 
winding by the d-q coordinates system respectively. We can write the relationship of electromagnetic torque 
as follows: 


Te = Pi (igs Yia ~ las Yiq) (3) 


where: P1 is the pair of poles for torque windings. 

In this work the field-oriented control system (FOC) is used of the BIM speed control, this system 
also enables the flux and torque of BIM to be controlled separately as in separately excited DC motor, in other 
words, the electric torque is controlled by the ig; current and the flux by the iias current. By using rotor field- 
oriented control scheme, the vector of rotor field for torque winding is aligned in the direct axis, that means 
the linkage flux components of the rotor for torque winding will be: Yar = 0 and Yar = Yr [2], [17]. 


* 
* Yir 


T a) 
a SO LırTe* 
tqs = PiLım Wir" (5) 
Depending of the machine parameters the slip speed can be write as (6). 
LimRiriisq” 
Ws, = Sam arasi (6) 


Lir Yir” 


From the measured speed of the rotor (@, and slip speed (@.1) can be obtained the flux transformation angle (0) 
as follows: 


0 = f (wsı + @,)dt (7) 


from basic equations of rotor. 
I 
Wy = lds (1 + tp) (8) 
Lm 


The electromagnetic torque equation can be determined as: 


L 
Te =P = (igs Wr) (9) 
r 


Li 


where: Liris the inductance of the rotor for torque winding, Riris the resistance of the rotor for torque winding, 
p is the partial derivative, t is rotor time constant, TL is the load torque, J and P; are moment of inertia and 
pole pairs for torque winding respectively. 


2.2. Suspension forces of bearingless induction motor 

Two-poles suspension force winding and four-pole torque winding are embedded in stator slots of the 
BIM. The electromagnetic coupling of the bearingless induction motors is very complex, because the coupling 
between the two-pole winding and four-pole windings, and the couplings between the stator windings 


Design of suspension control system for bearingless induction motor using ... (Qasim Kadhim Jasim) 


774 m) ISSN:2088-8694 


themselves. So, to deal with this form easily, the three phase (abc) is transformed to a two-phase system (aß) 
and then transformed to a two-phase rotation system (dq). Let W2a, W», Wia and W are the flux linkage of 
windings: 2a, 2b, la and 1b respectively and izas, i2qs, i1ds and igs are the stator current components in d-q 
coordinates of windings 2a, 2b, la and 1b respectively. The flux linkage equations of the BIM can be written 
as follows [18]-[21]: 


[¥] = [L] + [i] (10) 
Pia Lia Miaib Miaza Miazb has 
Pip z Minia Lip Minza Mibzb , | 9 (11) 
Poa Mzaa Mza1b Laa Mza2b lads 
Pop Mbia Mzp1b Mzp2a Lap l2qs 


where: M and L is the mutual and self-inductance and subscript of 1 and 2 refer to torque winding and 
suspension windings respectively. The value of mutual inductance between four- pole windings with 
themselves is zero, the value of mutual inductance between the two-pole windings with themselves is zero, the 
value of self-inductance for torque windings (Lis) and the suspension windings (L2;) are constant, only the 
value of mutual-inductance between the torque winding and the suspension windings is variable and dependent 
of the rotor displacement variation in X and Y direction, So the equation of inductance matrix for the (BIM) 
can be written as [20]: 


las 0 —MX MY 

0 Li, MY MX 
-MX MY L,, 0 

MY MX 0 Ly 


[L] = (12) 


where M is mutual inductance coefficient for fore-pole windings and two-pole windings, X and Y are radial 
displacement of rotor in X and Y axis respectively. According to the principle of energy conversion in the BIM, 
we can describe the magnetic energy that is stored in the stator windings as follows [21]: 


[w] = $+ [ATI (13) 


By neglecting magnetic flux saturation, the axial forces Fx = dw/dx in X direction and Fy = dw/dy in 
Y direction as follows: 


dw . 2 A 
F =~ .f-i i i 
p 4 l = m| ds | Fa (14) 
y dy. ligs lids 2qs 
Let, iias = Im * cos(wt) (15) 
igs = Im * Sin(wt) (16) 


then 


Fx] Hia —cos(wt) sin(wt)] lizas“ 
pa ms m 


sin(wt) cos(wt)] |izqąs 
where: Im is the magnetizing current of the torque windings. 


2.3. The radial motion equations of the suspended rotor 

In this section, we will study the axial (radial) motion of the rotor based on the newton's second law, 
which states that the sum of the forces that act on the rotor are equal to the product of mass multiplied by 
acceleration, the radial motion equations of rotor shaft in the X and Y directions as shown in Figure 3 can be 
written as [1], [22]: 


Fdx + Fsx+Fx = mË (18) 


Fdy + Fsy + Fy = mY (19) 
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Fsx = Ks X and Fsy = Ks Y 


_ lraB? 


Ks = 228 (20) 


where: m, l and r are the mass, length and radius of rotor respectively, Fsx, Fsy are unilateral magnetic pull 
components along x and y direction respectively, Ks represents radial displacement stiffness factor, B is 


magnetic-flux density, g is the length of air-gap, ois the vacuum permeability and Fdx, Fdy are the external 
disturbance of the rotor in X and Y axis direction respectively. 


(b) 


Figure 3. State space of rotor motion in (a) X axis and (b) Y axis 


3. THE PROPOSED CONTROLLER (FUZZY-D 

The proposed control system is simulated by using the MATLAB-Simulink platform as shown in 
Figure 4. This control system is divided into two subsystems for rotation torque control and suspension control. 
The rotation of the motor is controlled by using the field-oriented control system (FOC), while the axial 
suspension forces in the X and Y axis are controlled by using the differences between the references and 
measured values of rotor position, which are adjusted by using the fuzzy-I controller. The adjustable suspension 


currents izq* and izg* are obtained by converting the forces Fx and Fy to currents using the equation of 
modulation as (21) [20]. 
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Figure 4. Block diagram of the presented control system 


Design of suspension control system for bearingless induction motor using ... (Qasim Kadhim Jasim) 


7716 o ISSN:2088-8694 


By using dq-abc transformation the three reference currents of suspension force are determined. The 
suspension windings are supplied by three phase bridge inverter and the current of each phase is controlled 
separately, this is simply achieved by comparison the reference currents (is¥) to the actual currents (is) by using 
the hysteresis comparator as shown in Figure 5. In the bearingless induction motor, a new set of windings has 
been added with the main stator windings. The purpose of these additional windings is to keep the rotor in the 
center and away from the stator, which are controlled by the separate control system. 


SIX 
IGBT 


INVERTER 


Figure 5. Current controller of induction motor 


Many well-known structures are employed in design of the PID-like fuzzy logic controller depending 
on the structure used, fuzzy-I control system has been proposed in this paper. The structure is fuzzy logic 
connected in parallel with integration. This controller combines the features of both fuzzy and integration. 
Figure 6 illustrates the diagram of the fuzzy-I control system; two of these controllers are used, one for 
controlling the X position of the rotor and the other for Y position. 


K1 
g n A o > 
Diab et ply E 
K2 
Figure 6. The simulation of fuzzy-I controller 


K1 and K2 are selected based on the limits of error and the change in error to be suitable to the range 
of the fuzzy logic input, K3 and K4 are selected to improve the control system performance by using trial and 
error. The main structure of the fuzzy logic system includes three stages, the fuzzification, Inference rules, and 
defuzzification stage. In fuzzification, the data of input are converted to the fuzzy values, which are lingual 
values understandable by inference mechanism, where fuzzification is the process of applying membership 
functions to input variables. The inference rule employs fuzzy operations to map the input linguistic value to 
the output linguistic values. In the defuzzification, the fuzzy values which are resulting from the inference are 
transformed to crisp values [23], [24]. 

The fuzzy logic employed in this study consists of two inputs, the error of displacement (E) and change 
in error (CE), each one of the inputs and outputs has a membership function. The input and lingual variables are 
NB (negative big), NM (negative medium), NS (negative small), Z (zero), PS (positive small), PM (positive 
medium), and PB (positive big) where, the universe(rang) of both error and the change in error is [-1 1] as shown 
in Figure 7. There are six different kinds of membership function which are trapezoidal, triangular, Gaussian, 
sigmoidal, s-type and z-type. For input: triangle, s-type and z-type are used, while for output, triangle shape is 
used, this combination gives most accurate results compared with other kinds of the membership functions. 

The lingual variables of the output are NB (negative big), NM (negative medium), NS (negative 
small), Z (zero), PS (positive small), PM (positive medium), and PB (positive big) and their input 
range(universe) of [-1 1]. The arrangement of the output membership functions is shown in Figure 8. The 
fuzzy rules are statements set of IF-THEN used to determine the output state from the input state [25], the rules 
of the fuzzy control system are tabulated in Table 1. 
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Figure 7. Input membership functions (a) E (error) and (b) CE (change in error) 
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4. RESULTS OF SIMULATION 

To effectively verify the proposed controller, a system of bearingless three-phase induction motor is 
simulated and analyzed at different operating conditions. In this section the results of the proposed control system 
are compared with conventional PID results. The rotor deviation should not exceed 50 percent from the length of 
air-gap (0.6 mm) to prevent the collision with the stator. The parameters of BIM are tabulated in Table 2. 


Table 2. BIM parameters [26] 


Symbols Value Description Symbols Value Description 
Lm1 0.1586 Mutual inductance for torque-windings (H) m 2.86 Mass of the rotor (Kg) 
Lir 0.0092 Leakage inductance of the rotor (H) N1 60 No. of turn for torque windings 
Rir 11.48 Resistance of rotor (Q) N2 140 No. of turn for suspension windings 
J 0.00796 Moment of inertia (Kg.m°) P1 2 Pole-pairs for torque windings 
l 105 Long of the rotor (mm) P2 1 Pole-pairs for suspension windings 
r 97.8 diameter of the rotor (mm) g 0.6 Air-gab (mm) 


The motor is started at reference speed of 3000 rpm and load torque TL=2 Nm and no external 
disturbance forces on rotor shaft. At time 0.2 s, an external-disturbance force in negative X axis direction of 40 N 
is applied on the rotor as shown in Figure 9(a). Figure 9(b) shows the rotor deviation from the center point in the 
X axis direction at time 0.2 s. The value of the maximum rotor deviation is 160 um by using a conventional PID 
controller and 88 um by using proposed fuzzy-I controller and the peak-to-peak rotor deviation is 160 um by 
using a conventional PID controller and 123 um by using the proposed fuzzy-I controller. at time 0.2 s the radial 
suspension force (Fx) will increase from zero to 40N to overcome the external disturbance force as shown in 
Figure 9(c). Figure 9(d) shows the peak to peak three phase current of suspension windings, before the disturbance 
force is applied, the suspension windings current is 90 mA to counter the rotor weight and then at time 0.2 s the 
suspension current is increased to 145 mA to counter both of the rotor weight and the external-disturbance force. 

To verify the suspension control system performance during the variation of reference speed, at time 
0.5 s the reference-speed (@*) is increased from 3000 to 4000 rpm (314 to 419 rad/s) as shown in Figure 10(a). 
The speed of the rotor (wr) reaches a steady state after 27 ms. Figure 10(b) shows the effect of the speed 
variation of rotor position in X axis direction, the value of the maximum rotor deviation is 100 um by using a 
conventional PID controller and 53 um by using proposed fuzzy-I controller and the peak-to-peak rotor 
deviation is 138 um by using a conventional PID controller and 95 um by using the proposed fuzzy-I controller. 
Figure 10(c) shows the effect of the speed variation of rotor position in Y axis direction, the maximum value 
of rotor deviation is 78 um with traditional PID, this deviation decreases to 40 um by using the proposed 
fuzzy-I controller. Figure 10(d) shows the effect of speed variation on the suspension current, the motor current 
increases during the acceleration period, which leads to a decrease in suspension current from 145 mA to 
30 mA (peak to peak value) because the suspension current is inversely proportional to the load current in the 
event that the suspension force is constant. When the speed reaches to the steady state, the suspension current 
returns to the normal value (145 mA). 
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Figure 9. Performance of control system during disturbance force: (a) disturbance force, (b) X- displacement, 
(c) X-suspension force, and (d) the suspension current 
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Figure 10. Effect of the speed change: (a) the rotor speed, (b) X displacement, (c) Y displacement, and 
(d) the suspension current 


In case of speed and torque variation the displacement value towards the X axis is greater than the 
displacement towards the Y axis, because the displacement towards the X axes is affected by the force applied 
towards the X axis (disturbance force) which is equal to 40N, and the displacement towards the Y axis is 
affected by the rotor weight (28.5 N). 

To verify the suspension control system performance during the variation of load torque, at time 
0.9 s the load torque is increased from 2 to 10 Nm as shown in Figure 11(a). Figure 11(b) shows the effect of 
load variation of rotor position in X axis direction, the maximum rotor deviation is 43 um by using a conventional 
PID controller and 30 um by using the proposed fuzzy-I controller and the peak-to-peak rotor deviation is 43 um 
by using a conventional PID controller and 39.1 um by using proposed fuzzy-I controller. Figure 11(c) shows 
the effect of load variation on rotor position in Y axis direction, the maximum rotor deviation is 32 um using 
PID, this deviation decreases to 22 um by using the fuzzy-I controller. Figure 11(d) shows effect of the load 
variation on the suspension current, the load current increases with the increase of load, which leads to a 
decrease in suspension current from 145 mA to 103 mA (peak to peak value), because the suspension current 
is inversely proportional to the load current when the suspension force is constant. 

Summary of results and a comparison between the conventional PID and the fuzzy-I controller based 
on the peak-to-peak rotor deviation in X axis direction are shown in the Table 3. The percentage deviation is 
calculated by (22). 
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peak to peak deviation with PID—peak to peak deviation with Fuzzy_I 


Percentage deviation % = (22) 


peak to peak deviation with PID 


The simulation results show that the peak-to-peak rotor deviation by using fuzzy-I controller is decreased 
compared with conventional PID so the fuzzy-I controller is more suitable to control the radial suspension 
system of bearingless induction motor. 
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Figure 11. Effect of load torque variation: (a) load torque, (b) X displacement, (c) Y displacement, 
and (d) the suspension current 


Table 3. Rotor deviation in X axis with effect of different conditions 


Disturbance force effect 


Speed variation effect 


Load torque variation effect 


Maximum deviation with PID 160 (um) 100 (um) 43 (um) 
Maximum deviation with fuzzy-I 88 (um) 53 (um) 30 (um) 
Peak to peak deviation with PID 160 (um) 138 (um) 43 (um) 
Peak to peak deviation with fuzzy-I 123 (um) 95 (um) 39.1 (um) 
Percentage decrease of deviation 23 31 9 


CONCLUSION 
The rotor of the BIM is supported by magnetic forces, the modelling of these forces is complex. In 


general, approximation and simplification are used to create a linearized model in order to assess the 
performance of the proposed controllers. The proposed suspension controller is based on fuzzy-I controller, 
the fuzzy logic controller is nonlinear system and doesn’t require exact mathematical models and the 
integration helps to reduce the steady state error. The simulation of the suspension control system shows that 
the fuzzy-I controller reduces the peak-to-peak rotor deviation by 23% under effect of external disturbance, 
31% under effect of speed variation, and 9% under effect of load variation, this percentage decrease of deviation 
for the rotor position were found through the peak-to-peak deviation values. The propitious results certify the 
applicability of fuzzy-I to control the radial suspension system of bearingless induction motor. 
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